APS/123-QED 
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We study spin-1 bosons in an optical lattice under a magnetic field by the Gutzwiller approx- 
imation. Our results thus obtained join the discontinuous phase boundary curves obtained by 
perturbative studies through a first-order transition. On the phase boundary curve, we also find a 
peculiar cusp structure originating from the degeneracy between different spin Mott states under a 
magnetic field. The magnetic field dependence of both fluctuation in the total number of bosons 
and the spin magnetization clarifies that the superfluid phase is divided into two regions reflecting 
the coexisting first- and second-order superfluid transitions. 
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Superfluid (SF) and Mott insulators (MI) are very in- 
teresting subjects in condensed matter physics. Recently, 
the transition between the SF and MI of spinless Bose 
atoms has been experimentally demonstrated in an opti- 
cal lattice system [l| closely following the theoretical pro- 
posal based on a Bose-Hubbard model (BHM) 0,0 due 
to Jaksch et al. Q . In contrast to the case of fermions [j| , 
the Gutzwiller approximation (GA) |(| describes an ac- 
curate second-order SF-MI transition of spinless bosons, 
which is supported by quantum Monte Carlo simulations 


Trapped spinor bosons have also been investigated 
both theoretical ly |3| a nd experimentally Q. In an op- 
tical lattice, MI|iriTll| and SF[l(J H2 phases of spinor 
bosons have been studied. MI phases under a magnetic 
field have also been studied recently 01 . The SF-MI 
transition in a spin-1 BHM has also been studied by 
a perturbative mean- field approximation (PMFA) [lij . 
which treats the hopping process between adjacent sites 



as a perturbation [Q, by the GA [T^, by a conven- 
tional mean-field approximation (CMFA) [lj|, and by 
density-matrix renormalization group in one dimension 
[l8|. The GA and CMFA results show a possible first- 
order transition (FOT) at a part of the phase boundary 
curve, where the critical value of hopping matrix element 
just on the phase boundary is smaller than that obtained 
by the PMFA, which describes a second-order transition 
(SOT). On the other hand, the results obtained by the 
PMFA can also be obtained by the GA, which assumes 
a partial set of the full wave functions [l(j. Hence, gen- 
erally speaking, the results obtained by the GA improve 
those obtained by the PMFA. 

The PMFA approach has also been applied to the SF- 
MI transition in the spin-1 BHM under a magnetic field 
[T^.l20| . However, the phase boundary curve obtained by 
the PMFA in the limit of a weak magnetic field (B — > 0) 
[H i3 does not agree with that obtained by initially 
assuming zero magnetic field (B = 0) when the MI 
phase has an odd number of bosons. Moreover, the phase 
boundary obtained by the PMFA [l^. |20| is a discontin- 



uous function of magnetic field. These features might be 
attributed to the PMFA itself. 

In this paper, we study the spin-1 BHM by the GA. 
In contrast to the PMFA, the GA shows that the SF-MI 
phase boundary is a continuous function of magnetic field 
even around a zero magnetic field. In the phase diagram, 
we also find a special point where a degeneracy of MI 
states with different spins plays an important role. We 
also investigate superfluid properties in terms of mag- 
netization and fluctuation in the total number of bosons 
(FTNB) in the system. Both quantities, which are exper- 
imentally observable, have interesting magnetic field de- 
pendence originating from the coexisting FOT and SOT 
in the present system. 

The BHM of spin-1 bosons [H El III 13 E HI is 
given by H = Hh op + H mt + H m3bg in standard notation 
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Here, for simplicity, we assume that a magnetic field is 
parallel to the z axis and that the system is uniform. In 
this paper, we use Uo as a unit of energy and assume 
an antiferromagnetic interaction U2 = 0.04, which corre- 
sponds to 23 Na atoms. 

The Gutzwiller variational wave function for a site is 
defined as $ = Z N g(N)\N), \N) = £ s f(N, S)\N,S), 

and \N,S) = 52sA N > S >. S *)\ N > S > S *)- Here ' \ N > S > S *) 
has N bosons, a total spin S, and a magnetic quantum 
number S z , where S must be odd (even) for an odd (even) 
N [TfJ. The variational parameters must satisfy the 
normalization condition ^2 N ^(AQI 2 = J2 S \f(N, S)\ 2 = 
J2s S, S z )\ 2 = 1. In this paper, we take the 

complete set from N = to N = 6, which is suffi- 
cient for a numerical convergence. We employ a stan- 
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FIG. 1: Phase boundary curves as a function of chemical 
potential u. SF and MI indicate the superfluid and the Mott- 
insulating phases, respectively. The solid and dashed curves 
are obtained for b = 0.001 using the GA and the PMFA, 
respectively. The dot-dashed curves around the MI phase 
with TV = 1, 3 are obtained by the PMFA under zero magnetic 
field and are not visible around the MI phase with N — 2,4 
because they are indistinguishably close to the dashed curves. 



dard definition such that the MI phase has zero FTNB 
((N tot - (N tot }) 2 )/(N to tUN tot = while the SF 

phase has a finite one |2lJ. In an SF phase close to an 
MI phase with N bosons, we consider g(M) for M ^ N 
as an SF-order parameter because finite g(M) results in 
finite FTNB. The SF order parameters have finite (no) 
jumps at the SF-MI phase boundary for a FOT (SOT). 

Superfluid-Mott-insulator transition — Figure 1 shows 
phase boundary curves on the zt-[i plane (z: the number 
of adjacent sites) under a very weak magnetic field b = 
0.001. In addition to the results obtained by the PMFA 
and the GA, we also plot those obtained by the PMFA 
under zero magnetic field (dot-dashed curves). 

We can easily see that the results obtained by the 
PMFA for b = and b = 0.001 are clearly different from 
each other around an MI phase with an odd number of 
bosons. This difference originates from the calculation 
procedure. Namely, a degenerate PMFA is employed to 
lift the degeneracy among S z = ±1,0 states under zero 
magnetic field , while the degene racy has already been 
lifted under a finite magnetic field [n| |2(| • Because the 
PMFA neglects low-spin states even under a weak mag- 
netic field, the PMFA overestimates the antiferromag- 
netic interaction energy in a possible SF phase, resulting 
in a large critical value of zt for the SF-MI transition. 
The results obtained by the GA fall in between the two 
results obtained by the PMFA. The SF-MI transition is 
a FOT at a part of the phase boundary, where the re- 
sults obtained by the GA do not completely agree with 
those obtained by the PMFA. In the limit of b — ► 0, the 
results obtained by the GA are not below those obtained 
by the PMFA initially assuming 6 = 0, but completely 
agree with them. This is consistent with the previous re- 
sult under zero magnetic field 0] such that a FOT only 



FIG. 2: Phase boundary curves on a b-zt plane for fi = 2.5 
and for /i = 1.85. The dashed (solid) curves are obtained by 
the PMFA (GA). The inset is an enlargement around b = 0.06 
for u = 1.85. 



occurs for very small U2 around MI phases with odd N 
bosons and does not occur for any U2 when N = 1. 

On the other hand, around the MI phases with an even 
number of bosons, the phase boundary curves obtained 
by the GA (the PMFA) for b = 0.001 are almost indis- 
ting uishably close to those for b — 0, which are from Ref. 
[lfj (Ref. [13 )• This shows that the singlet MI phases 
are robust under a weak magnetic field. 

Figure 2 shows the magnetic field dependence of the 
phase boundaries for \i = 2.5 and /i = 1.85. For /1 = 
2.5, the PMFA assumes an MI state with = 3 has a 
spin (S = 1 or S = 3) that discontinuously depends on 
whether b is smaller or larger than 0.1 j^. This results 
in a curious jump of the critical value of zt for the SF- 
MI transition just for b = 0.1. The solid curve obtained 
by the GA gives a continuous phase boundary by joining 
the separated dashed curves and exactly agrees with the 
dashed curves except for a finite region 0.0838 < b < 0.1, 
where the transition is a FOT. 

On the other hand, for [i — 1.85, the solid curve 
obtained by the GA has a sharp cusp structure around 
b = 0.06 where the MI states with S = and 5 = 2 are 
degenerated. The critical value of zt just for b = 0.06 
is different from both of the two limiting values from 
the weak or strong magnetic field region obtained by 
the PMFA. The transition is a FOT near (but not 
just for) b = 0.06. However, the SF order parameters 
become smaller when b becomes closer to 0.06, and 
finally the transition becomes a SOT just for b = 0.06. 
In fact, by using a degenerate PMFA that assumes 
two MI states such that \N,S - 2) and \N,S) as 
zcroth-order states, we obtain the critical value of zt 

/ , \ — l 1 a-c 1 1 / 1 1 a - \ / c 2 — 4b-b + 

as (zt deg ) - 5 |°+ _ 2b~ + 2V 1 + ]b~\)y a'i-b. J' 
where a± = (a Sj+ ± a s _ + a s _ 2 ,+ ± «s , -2,-)/2, 
b± = [as,+ ±«s - - (as-2,+ ±as-2,-)] 2 /4±/3 2 , and c = 
[(as,+ - as-2,+) 2 - (as,- ~ «S-2,-) 2 ]/2. Here 
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E(i,j,k) is the energy per site of an MI state 
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U i h\ and f l 'i' k — \{l,rn,n\~a\i,],k)\- 
\l,J,K) ana J lttni7l = E {i,j,k)-E{l,m,n) ' 

a creation or an annihilation operator which 
\l,m,n) and \i,j,k) [23| . The PMFA chooses 
as the critical value only when zt^ eg is smaller 
than the two limiting values from the weak or strong 
magnetic field. Furthermore, another condition 

[c - asgn(b-)yj((? -4b-b+)/(a 2 _ - fe_)]/(26_) <1 

must also be satisfied because the absolute values of the 
SF-order parameters must be non-negative. The latter 
condition is not satisfied in the case of /j, = 2.5 in Fig. 2. 
It should also be noted that the same critical value as 
zt^eg can be not only numerically but also analytically 
obtained by the GA including the states that emerge 
as zeroth-ordcr or intermediate states in the degenerate 
PMFA calculation. 

Superfluid properties — Figure 3 shows magnetization 
per site and FTNB as a function of magnetic field for 
fj, = 1.92 and for a constant zt = 0.02 24]. The FTNB is 
proportional to the experimental observables, such as the 
compressibility and the inverse square of the sound ve- 
locity [2| ■ Both the magnetization and the FTNB curves 
can be clearly divided into four parts depending on the 
magnetic field. There are always discontinuous jumps of 
the differential magnetic susceptibility and the derivative 
of the FTNB on the boundaries between the four parts. 

To clearly understand the magnetization and the 
FTNB curves, we also show phase boundary curves in 
Fig. 4, where the SF-MI transition is a SOT (FOT) un- 
der a magnetic field b < 0.06 or 0.0670 < b (0.06 < b < 
0.0670). The MI states \N = 2, S = 0} and \N = 2, S = 
2) are degenerated at b — 0.06 and the former (latter) is 
more favored under a weaker (stronger) magnetic field. 
For zt = 0.02 in Fig. 4, the first region (b < 0.0462) 
and the fourth region (0.0612 < 6) correspond to the MI 
phases \N = 2, S = Q) and \N = 2, S = 2), respectively, 
where the magnetization is constant and the FTNB is 
zero as shown in Fig. 3. 

The SF phase in the second region 0.0462 < b < 0.0596 
has a perturbative character. The SF state is continuously 
connected to the nearest MI state, and the spin property 
of the system is close to that of the MI state. The in- 
set of Fig. 4 indeed shows that |/(2,2)| 2 , which is the 
amplitude of a high spin state \N = 2,S = 2} normal- 
ized as |/(2, 0)| 2 + |/(2, 2) | 2 = 1, is almost negligible and 
that the lowest spin state \N = 2, S = 0} is dominant in 
the second region. This perturbative character originates 
from a SOT between the SF phase in the second region 
and the MI phase in the first region. 

On the other hand, the SF in the third region has 
a non-perturbative character such that states with high 
spins are largely included, resulting in the large (S z ) as 
shown in Fig. 3. The inset of Fig. 4 shows that |/(2, 2)| 2 
is indeed large in the third region. This change of the 
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FIG. 3: Magnetization per site (S z ) (solid curves) and fluctu- 
ation in the total number of bosons {{Ntot - {N tot )) 2 )/(N tot ) 
(dashed curves) as a function of magnetic field for fi = 1.92. 
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FIG. 4: Phase boundary curves for \i = 1.92. The solid 
(short dashed) curve is obtained by the GA (PMFA). The 
long dashed line, which shows the zt — 0.02 assumed in Fig. 
3, intersects with the solid curve at b = 0.0462 and b — 0.0612. 
The inset shows |/(2, 2)| 2 as a function of magnetic field (See 
text for the definition). 



SF character also affects the FTNB as shown in Fig. 3. 
It should also be noted the non-perturbative SF phase 
can also be characterized as a coherent- state-like charac- 
ter with a large kinetic energy as in the case of a zero 
magnetic field [ifil ]. This non-perturbative character is 
related to a FOT between the SF phase in the third re- 
gion and the MI phase in the fourth region. 

We can see the crossover between the two SF phases 
in a wide parameter region when the SOT and FOT 
phase boundary curves coexist. For instance, the critical 
value of zt for the SF-MI transition just at the bound- 
ary between MI phases with different spins has not to 
be necessary a local minimum as a function of magnetic 
field. Although not shown here, we indeed found a clear 
crossover between the non-perturbative and perturbative 
SF phases for zt > 0.075 and fx — 2.5 around b = 0.1 (See 
Fig. 2 for the phase boundary curve). Here, the critical 
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FIG. 5: The same plot as in Fig. 3 for fj, = 1.85 and zt = 0.07. 



value of zt = 0.075 for b — 0.1 is not a local minimum as 
a function of magnetic field. 

Let us explain some details. There are no gaps on 
the curves of the magnetization and the FTNB at the 
boundary between the perturbative and non-perturbative 
SF phases in Fig. 3. On the other hand, there are gaps 



in the parameter region as shown in Fig. 5 for /i = 1.85 
and zt — 0.07, for which the phase boundary curves are 
shown in Fig 2. For (j, = 1.85 and zt = 0.07, the SF- 
MI transition occurs only once for b ~ 0.0274 through a 
SOT. However, a FOT also occurs for a slightly larger b 
and a slightly smaller zt, resulting in a crossover between 
the two kinds of SF phases for zt = 0.07 p5j . 

We finally note that the GA neglects inter-site cor- 
relation effects, which could be important when U2 is 
smaller than or comparable with zt 2 /Uo and/or when 
the dimension of the lattice is low. Although zt 2 /Uo is 
somewhat smaller or much smaller than U2 in the typi- 
cal parameter sets we have assumed in this paper 26], a 
somewhat larger U2 would more clearly justify the GA. A 
larger {/ 2 is also favorable from the experimental point of 
view. This is because both Zeeman energy and antiferro- 
magnetic interaction energy become comparable under a 
stronger magnetic field and the interesting features clar- 
ified in this paper can be easily observed. 
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